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Fermion family recurrences in the Dyson-Schwinger formalism 
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Abstract. We study the multiple solutions of the truncated propagator Dyson-Schwinger equation for a 
simple fermion theory with Yukawa coupling to a scalar field. Upon increasing the coupling constant g, 
other parameters being fixed, more than one non-perturbative solution breaking chiral symmetry becomes 
possible and we find these numerically. These "recurrences" appear as a mechanism to generate different 
fermion generations as quanta of the same fundamental field in an interacting field theory, without assuming 
any composite structure. The number of recurrences or flavors is reduced to a question about the value of 
the Yukawa coupling, and has no special profound significance in the Standard Model. The resulting mass 
function can have one or more nodes and the measurement that potentially detects them can be thought 
of as a collider-based test of the virtual dispersion relation E — y/p 2 + M(p 2 ) 2 for the charged lepton 
member of each family. This requires three independent measurements of the charged lepton's energy, 
three-momentum and off-shellness. We illustrate how this can be achieved for the (more difficult) case of 
the tau lepton. 

PACS. 12.15.Ff Quark and lepton masses and mixing - f f.30.Rd Chiral symmetries 
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1 Introduction 

Why Nature has laid down exactly three fermion families 
in the same representation of the Standard Model's gauge 
groups with universal couplings to the gauge bosons re- 
mains to be explained. Many ideas have been presented 
in the literature and we quote some to exemplify the far- 
reaching implications of any experimental progress in the 
first question of what has been called the "Fermion prob- 
lem". 

"Democratic approaches" [T] are in general based on 
the idea of equal fermion-to-Higgs Yukawa couplings. One 
then needs mechanisms to generate specific lepton or quark 
mass patterns, see for example [2J. Let us also recall the 
classic work of Froggat and Nielsen [3] that introduces the 
concept of " horizontal flavor" symmetry (nowadays often 
called Interfamily symmetry) by which the three families 
are degenerate or quasi-degenerate at some high scale be- 
yond present experimental reach, due to unknown symme- 
try. At lower energies dynamical (renormalization group) 
effects amplify any small symmetry breaking term causing 
the ratios between masses that we see in current acceler- 
ators. In essence, the idea is that each of the elementary 
particle types, for example the charged leptons (e,fi,r) 
provide a representation of the posited interfamily sym- 
metry. Since the symmetry is broken at the TeV scale 
and below, revealing it requires theoretical extrapolation 
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from current experimental data. The same authors have 
followed-up with ideas based on Anti-grand unification [4] . 

An alternative line of research assumes that the known 
elementary particles are really composite objects of more 
fundamental "preons" [5]. To date no experiment has re- 
vealed composite structure beyond the usual layers of quan- 
tum pairs involved in radiative corrections around seem- 
ingly point-like sources. One can argue that a new strong 
interaction binds preons and this might conceivably be 
revealed in future experimental efforts. 

Further, dynamical breaking of chiral symmetry as a 
genuine quantum effect is well known in strong interaction 
physics [51(7] and has also been invoked in flavor physics 
as an alternative to the tree-level Higgs mechanism, for 
example by Brauner and Hosek [5]. 

However no one seems to have paid attention to the 
fact that the Dyson-Schwinger equations do provide one 
with several solutions depending on the strength of the 
coupling. In this work we call attention to this point and 
its potential relevance for the fermion problem. Each of 
the three fermions corresponds in this hypothesis to a 
quantum over a different vacuum (but in the covariant 
formalism one pursues the study of correlation functions, 
here propagators, and sidesteps the issue of the vacuum 
wave functional). Far from attempting a complete theory, 
we will show the general features of the mechanism within 
a simple model of a fermion field coupled to a scalar boson 
through a Yukawa coupling. 
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From this prospective, different flavored fermions, one 
corresponding to each family, are quanta of the same field, 
so the lagrangian can be more economically written. Each 
fermion is an elementary one-particle excitation on top 
of a vacuum that is a local minimum of the Hamiltonian. 
The excitations over the ground-state vacuum provide the 
lightest fermion family. Finally, at fixed coupling g, the 
Dyson-Schwinger equation has a finite number of solu- 
tions, so a finite number of families arise. The question on 
why three families becomes only a question of the value 
of g, the interaction coupling, and may have no special 
significance. 

The equivalent of the propagator Dyson-Schwinger equa- 
tion in a non-covariant framework is the well studied mass 
gap equation of potential models of QCD. The equation 
has been numerically solved for the harmonic oscillator 
potential [9], [10], the linear potential [11] . and the lin- 
ear plus Coulomb potential [I21I11I- Bicudo, Ribeiro and 
Nefediev have systematically studied the excited solutions 
of this equation [15j . Our results seem compatible with 
this prior work. 



2 Dyson-Schwinger equation and its 
numerical solution 

We examine a simple Yukawa theory for one fermion field 
coupled to a real scalar field with Lagrangian density in 
Euclidean space 

C = W{i t>- m w )V + gm<j) + - m 2 )cj) ■ (1) 

The Dyson-Schwinger equation for the fermion propa- 
gator in this theory is represented in figure [T] in the rain- 
bow approximation (neglecting the vertex dressing). 

To examine generic features of dynamical chiral sym- 
metry breaking we can further ignore the running of the 
wavefunction renormalization and set Z(p 2 ) = 1, a con- 
stant. This leaves one scalar equation for the fermion mass 
function that reads 



M(p 2 ) = mvZ 2 



dqF(q,p)- 



M{q 2 



M 2 (q 2 ) 



In the chiral limit mq, — > this equation becomes homoge- 
neous and accepts a chiral-symmetry preserving solution 
M = that continuously deforms into a soft-running mass 
form when is not zero. But for strong enough kernels F 
it also admits other solutions that break chiral symmetry. 
The kernel in this simple Yukawa theory is 



F(q,p) 



47T 3 



\J\ - y 2 dy 



1 



{q-p) 2 + ' 



(3) 



with y the cosine of the polar angle in Euclidean four- 
dimensional space from p to q. 

To obtain excited solutions for M as function of p 2 
we apply an iterative linear method with different initial 
guesses with or without nodes for M(p 2 ). This method 



proceeds by examining linear deviations of the exact so- 
lution 

M(p 2 ) = M (p 2 ) + e(p 2 ) (4) 

that yield 



(5) 



m 9 Z 2 - M (p 2 )+ / dqF(q,p) 



M (q 2 



q 2 + M 2 (q 2 ) 



This can in turn be discretized (simply discretizing p) and 
cast as a linear system 



A tJ e 3 = (RHS)i 



(6) 



and solved with standard numerical linear algebra tools. 
The kernels require only the calculation of 2-dimensional 
integrals that are also standard fare in modern computers. 

By naive dimensional analysis one sees that the inte- 
gral needs UV regularization. We accomplish this by a 
simple cut-off method. Then the parameters should be 
chosen to run with the cutoff g(A), m^,(A) to ensure that 
the mass functions turn largely independent of A. In our 
results we display the dependence of the solution with g 
at fixed A. The boson mass is chosen to be 1 and sets 
the unit of the theory. 



3 Numerical results and discussion 

In this section we present the chiral-symmetry breaking 
conventional and new excited state solutions. By running 
our computer program incrementing g we meet critical g 
values go, g±, ...g n above which there are exactly n + 1 
solutions, the last with n nodes, that we find numerically. 
For example the first solution appears for a critical go — 
6.8 (of course this number runs with the cutoff). 

We show in figure [5] the three solutions obtained with 
g = U,A= 200. 

Then in figure [3] we plot the nodeless solution at fixed 
cutoff incrementing g sequentially to show the dependence 
of the mass function on the coupling. 

Only the nodeless solution has been widely used in 
past literature for its applications in hadron physics. The 
others are sometimes rejected [16] on the basis of "wrong" 
UV asymptotic behavior, arguing that they do not match 
with conventional mass running in perturbation theory, 
and artificially imposing a boundary condition that dis- 
cards all but one solution. 

We point out here that this mismatch poses no prob- 
lem. One should interpret these solutions as the one-particle 
excitations above different extrema of the Hamiltonian or 
vacuum replicae. Therefore they cannot all be matched to 
the same perturbation theory around one particular vac- 
uum. 

If one insists on viewing them from the one vacuum 
that connects smoothly to perturbation theory, then they 
have to be written down as complicated collective states 
in terms of the one-particle solutions over this vacuum. 
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Fig. 1. Dyson-Schwinger equation for the fermion propagator in Yukawa theory, in the rainbow approximation. 

Comparison of the three solutions 

g=14,m=l, A=200 




50 100 150 200 

p (arbitrary units) 



Fig. 2. We plot together the three solutions to the Dyson-Schwinger equation for a fermion-scalar Yukawa coupling g — 14, 
cutoff A — 200, and the scale set by the scalar mass being fixed to = 1. Note the equation for M in the chiral limit (zero 
current fermion mass) cannot distinguish M from —M so the global sign of any one solution in this graph is irrelevant and for 
display purposes only. In particular all M(0) can be taken as positive. 

In the particular case of the Yukawa theory, the an- 
gular integral in eq. [3] can be performed analytically |14j , 
but we prefer keeping a two-dimensional integral thinking 
of future work with more complicated model Lagrangians. 

In agreement with the findings of [15], we see that the 
higher excited vacua in this covariant Euclidean formu- 
lation present zeroes in the mass function. A version of 
the Sturm-Liouville theorem must be at work (this is not 
surprising since the Hamiltonian is hermitian). 

Also note that the values of M(0) generated in this 
covariant Yukawa theory show quite some hierarchy. In 
figure [2] one can see that they are well-spaced. Therefore 



this field-theory mechanism might conceivably be at work 
in the fermion flavor problem. However M(0) seems to be 
larger for the nodeless solution than for the excited ones 
(this puzzling property also appears in the equal-time ap- 
proach of [15], and the extent of its model dependence 
needs to be further investigated). To confirm this point, we 
choose a different set of parameters and implementation 
of the computer code. With = 100, a cutoff A = 10 4 
implemented as a Gaussian fall-off of the integration mea- 
sure (as opposed to terminating the grid), and g = 50, 
with a small bare fermion mass = 5 • 10~ 4 , one ob- 
tains the three solutions depicted in figure [4] This graph 
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Fig. 3. Nodeless solution for increasing values of the coupling constant (m^> = 1 and A — 200 are held fixed). Note how M(0) 
grows steadily with g (the sign of M is irrelevant). 



confirms that the nodeless solution is higher in mass, and 
that it is possible to generate sizeable mass hierarchies 
between the solutions. We use this run to also illustrate 
the dependence of M 2 with the cutoff, displayed in figure 
[H As can be seen, the excited solution does not decrease 
monotonically as g — > 0, it rather seems to become a sin- 
gular solution at the critical coupling. This technicality is 
however probably irrelevant for physical applications. 

In hadron physics the strength of the interaction being 
fixed by extensive phenomenology, there would be no free- 
dom to vary g. In the physically relevant cases studied in 
[T5] , excited solutions exist at physical values of the cou- 
pling. In [T7] Ribeiro and Nefediev have further extended 
their original results to include bound states (mesons) con- 
structed with the replicated quasiparticles. In particular 
there are Goldstone bosons over the replicated vacuum, 
that are of course not massless pions but would appear as 
excited pion states (separated from the ground state pion 
by the mass gap between the two vacua). Whether any 
conventional mesons accept a more convenient description 



as replicated mesons instead of excited mesons over the 
standard vacuum remains an open question. 

Working in Euclidean space as we do provides fastly 
convergent integrals. This technicality may be avoided by 
employing a Lehmann representation |18j that allows a 
direct solution in Minkowski space with similar results. 

4 Experimental signature 

In a free field theory, a field quantum of mass m satisfies 
the equation 

E = VP 2 + "i 2 (7) 
that causes a pole in the free propagator 



In an interacting field theory this pole will be shifted 
due to renormalization from its bare to its physical po- 
sition. The dependence on the renormalization scale can 
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Multiple gap functions solving the same DSE 

(shifted for clarity) 
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Fig. 4. Three solutions to the gap equation in Yukawa theory obtained with m^> = 100, A — 10 4 , = 5 ■ 10~ 4 , g = 50, 
confirming that the nodeless solution is higher in mass, and that a sizeable hierarchy of masses is possible. To represent them 
in a log-log plot, Mi and M2 being negative in part of their domain, the solutions have been shifted by a constant amount to 
positive values. 



be traded by a dependence on the particle Euclidean 4- 
momcntum so that 



S 



(9) 



p 2 - M(j> 2 ) 2 

yielding the transcendental equation for the pole position 



E = Vp 2 + M{E 2 -p 2 ) 5 



(10) 



Another way of visualizing the function M(p 2 ) of the 
interacting theory is to trade p 2 by the space-like part 
p 2 as argument of M, generating therefore a non-trivial 
dispersion relation for virtual particles, and the zeroes of 
the mass function M appear as points where E = |p| that 
would otherwise only be reached asymptotically at large p 
(see figure [7]). Tests for a non-trivial (real particle) disper- 
sion relation have been proposed and some carried out in 
the search for violations of Lorentz invariance [TH] • How- 
ever these are very indirect and usually performed at low 



energy, so an accelerator-based test of the (virtual parti- 
cle) dispersion relation by separately measuring E, p and 
the off-shellness is preferable. In our case, real particles 
still provide representations of the Poincare group and 
have constant, physical mass. However virtual particles 
off their mass shell will display the running mass function 
and this can be captured by analyzing the amplitudes of 
physical processes in perturbation theory. 

The corresponding dispersion relations are plotted in 
figures [H and [7] 

Lacking a real theory there are two possibilities we can- 
not discern at present. Either all three families of fermions 
arise as solutions spontaneously breaking chiral symme- 
try, in which case two of them should already have a zero 
in their mass function, or the first generation reflects the 
trivial solution only slightly modifying the bare mass, and 
the second and third generation solutions with dynam- 
ically broken chiral symmetry, in which case only one of 
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Dependence of M with the coupling g 
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Fig. 5. Dependence of M2 with the coupling constant at fixed m^> = 100, A = 10 4 , = 5 • 10 -4 . Upon comparing with figure 
[3]we see that the qualitative behavior when the solution approaches the critical coupling and is about to cease existing, here 
g — 43, is quite different for the ground and excited solutions. 



them would present a node. In cither case a zero is present 
in the mass function of one of the fcrmions. 

An interesting experiment is therefore a measurement 
of p, E and the off-shellness A 2 = p 2 — M 2 {p 2 ) of a r, 
at various momenta p. The B factories have accumulated 
large r samples [23] that are under analysis and can be 
employed for this purpose. 

The process one may investigate is depicted in Fig. [10] 
in the appendix below. An electron-positron pair collides 
and annihilates into a t~ t + pair, where the r~ is off its 
mass shell. Its four-momentum can be infered from the 
center-of- momentum energy of the incoming e~e + pair 
and by reconstructing the energy and momentum of the 
on-shell r + from its decay products, as detailed shortly. Fi- 
nally, one needs the off-shellness, but this can be obtained 
from the number of counts with given E, |p|. This is ob- 
vious since the cross section for a specific process where 
the off-shell t~ decays into for example l~viv T , depends 
on the off-shellness A 2 of the r~ (see also App. [XJ. 



Let us now detail the possible analysis: 

In the center of mass frame of an e~e + collision one 
can identify two back-to-back hadron jets, tagging the 
flavor by demanding that one of them kinematically 
reconstructs an on-shell r. This is rendered difficult 
by the undetected neutrino, that we assume to be the 
only unreconstructed track on the left side of the event. 
Complete reconstruction of the energy and momentum 
on the left side of the event (see Fig. [8]) is possible with 
a vertex detector (giving the direction of motion of 
the on-shell r) . The total energy is taken from collider 
calibration, and matched to the energy of the visible 
tracks. Balancing energy provides the energy of the 
missing neutrino, and automatically the momentum of 
the left side, that can be tested against the hypothesis 
of a physical r being produced. 

Next one examines the right side of the reaction, where 
the secondary and primary vertex seem to coincide (as 
appropriate for a virtual, off-shell r decaying rapidly). 



Llanes-Estrada, Van Cauteren and Paramo: Fermion family recurrences in DSE's 



7 



Dispersion relation for a fermion with running mass 

Arbitrary units 



2 2 1/2 

E(p)=(p +m ) 

E(p)=(p 2 +M(p) 2 ) 1/2 

E(p)=p (free massless theory) 
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Fig. 6. The virtual particle dispersion relation at fixed cutoff scale for the nodeless M(p 2 ) solution. Note that the interacting 
theory falls below the usual Yukawa dispersion relation because the mass function drops with p. This is not very different from 
standard perturbative running with p except the effect is much larger in the case of dynamical chiral symmetry breaking. 



The total three-momentum is known a priori. One bal- 
ances momentum to obtain the momentum taken by 
the undetected v T , and obtains the total energy. Now, 
it is not necessarily true that E 2 — p 2 = m 2 . 
4. The number of counts with given E, p on the right 
side is normalized to a cross section in terms of the 
off-shellness <i(A 2 ) and used to obtain M from the 
value of the intermediate propagator in perturbation 
theory. 

An alternative possibility is to carry out the measure- 
ment on the expected 0(20-30) r events at the OPERA 
experiment [21] in Gran Sasso National Laboratory. Given 
the momentum of the neutrinos from the CNGS beam 
one can attempt full kinematical reconstruction of (E, p). 
However here the difficulty resides in the off-shellness since 
a small value will be forced for the identification of the r. 

The experimental test for the muon and the electron is 
even simpler, for example through Compton and reverse 
Compton scattering where the fermion is off-shell in the 



intermediate state. For quarks one needs to take into ac- 
count that they always appear in bound states, and there- 
fore the mass function is always under an integral sign and 
reconstruction from experimental data difficult. However 
a particularly simple case is 

e~e + — > 7* — > ss 

where the strangeness be completely tagged by counting 
all hyperons in both jets. 

In principle the same analysis carried out here could 
also be undertaken for Majorana particles but this is out 
of our scope in this preliminary work. 



5 Conclusions 

We have called attention to an interesting feature of field 
theory, namely the possibility of having several solutions 
of the one-particle Dyson-Schwinger equations for a broad 
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Dispersion relation for a fermion with running mass 

Arbitrary units 




P 



Fig. 7. The virtual dispersion relation at fixed renormalization scale for the solution with two nodes (we find similar result 
for one node). As the running mass actually has a zero, the dispersive curve touches the massless limit that would otherwise be 
reached only asymptotically. 




Fig. 8. In center of mass frame production where p T + pr = one can fully reconstruct the decay of both r leptons if only 
the neutrinos escape undetected. The analysis steps are described in the text. 
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class of theories. The excited solutions appear upon in- 
creasing the coupling constant sufficiently, breaking chiral 
symmetry. 

The spectrum of states on top of each of the vacua 
is called "replica" or "recurrence" and provides a mech- 
anism that might be interesting for the theory of flavor. 
The number of recurrences arises from no special symme- 
try, since non-linear equations can perfectly have a finite 
number of solutions (again depending on the value of the 
coupling, this can be chosen to be three). 

These solutions have been shown in the past in a Hamil- 
tonian framework [15j and we have minimally extended 
the results of these authors to show that the covariant for- 
mulation allows similar phenomena. The excited solutions 
have mass functions that present zeroes and this reflects 
in the virtual dispersion relation (feature of an interacting 
field theory that should not be confused with breaking of 
Lorentz invariance). 

Whether these solutions will be relevant to the fermion 
family problem can be tested at B factories employing 
their large r samples directly produced. 

At this point it is a meaningless exercise to attempt 
to obtain the fermion masses in this scheme. Three data 
points Mj(0) should be used to fit three parameters, the 
current fermion mass (that we have here set to zero), the 
boson mass, and g at a given cutoff, so the current pre- 
dictive power is null. However we find this field theory 
mechanism still worth attention given the ease with which 
fermion coupling universality can be incorporated by hav- 
ing three recurrences of the same field. 

Testing the off-shell dispersion relation should allow 
to immediately discard the scenario so we believe we have 
presented a valid physical hypothesis. Of course, the suc- 
cess of radiative corrections in the standard model make 
one wish for a measurement at the highest possible energy, 
maybe at a future linear collider. 

This work has been performed in the framework of the 
research projects FPA 2004-02602, 2005-02327, PR27/05- 
13955-BSCH (Spain) and is part of the Masters thesis of 
Mr. Paramo Martm presented to the faculty of U. Com- 
plutense). TVC is a postdoctoral fellow for the Fund for 
Scientific Research - Flanders and acknowledges the sup- 
port of the "Programa de Investigadores Extranjeros en la 
UCM - Grupo Santander". 



A r-pair production in e + e collisions 
A.l On-shell cross section 

We look at the e + e~ — > t + t~ process depicted in Fig. [9l 
where all ingoing and outgoing particles are on their mass 
shell. Using the conventions of Ref. [22], the amplitude for 
the depicted process is 



iM=u T n -(p ) (-ieY) v T n t{p + ) 



Q 2 



x vU(k+) (tef) uf(k~) . (11) 



For the unpolarized cross section, one needs the squared 
amplitude averaged over initial and summed over final 
spins 



all spins 



x vU (fc + ) TV*- u f ik~)u%- (*-)7V«'+ 4+ ( fc+ ) 



x U n>-(P~) > ( 12 ) 

where the spinor indices are explicitly shown. For e + e~ 
annihilation, one has Q 2 = s. One can now rearrange the 
factors in Eq. (fT2)) . Making use of the relations 

u e(p)u e (p) = (j)+mi) ££ , , 



(13a) 

*4(p)«£'(p) = (J>- mi)#> > ( 13b ) 
one can write Eq. (fl"2|) as a product of traces 

1 e 4 



£|a<I s 



4 s 



2 Tr[(r 



1" {P 



m e )T 



x Tr [(J> + m T ) 7 P (f> + - m T ) 7„ 
Neglecting the electron mass, this results in 



(14) 



Ik ■ p k ■ p 



■ p k ■ p 



+ m 2 T k~- k + ) 

This can be written in terms of Mandelstam variables 

s = (AT - k+) 2 - (p+ - p-) 2 = Q 2 , 
t=(k- -p-) 2 = {p+ -k+) 2 , 
u= (AT -p+Y = (p- -k+) 2 , 



(15) 



(16a) 
(16b) 
(16c) 



Ei^i 2 = 



2e 4 



(* 2 



4 m 2 s 



2mi 



(17) 



where at high energies, the mass terms become negligible. 
In the center-of-momentum (com) frame, the amplitude is 



J2\M\ 2 = ^ [s 2 (1 + cos 2 (e T )) 

+ 4 m 2 s (1 -cos 2 (6» r ))] , (18) 



with T the corn-angle between incoming electron e~ and 
outgoing t~ . 

The differential cross section in the corn-frame is then 
given by 



da_ 



a 

4s» 



4 TO 2 



x [s(l + cos 2 (6» T )) +4m 2 .sin 2 (6g] , (19) 



with a = e 2 /47r. The total cross section is 



47ra 2 
3s 



Ami 



2mt 



(20) 
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Fig. 9. On-shell r-pair production in e + e collisions. The particles' momenta and spins are indicated, as well as the conventions 
for the photon propagator and the vertices. 




conventions for the boson propagators and the vertices. 



A. 2 Half off-shell cross section 

In this part of the Appendix, we look at the e~ e + — > t~t + 
process, where one of the r-particles is off its mass shell 
and decays into a lighter lcpton I (electron or muon). We 
discuss the case for an off-shell r~ . The case for an off- 
shell t + is clearly completely analogous. The process is 
displayed in Fig. flTJl 

The off-shellness of one of the r's changes the kinemat- 
ics. The total energy in the reaction is unevenly divided: 
E T + 7^ E T - ^ Wcom/2, where W COYa — yfs is the total 
energy in the corn-frame. The amplitude for the process 
depicted in Fig. [TO] is 

- iM = C (pi) (-i 9 -^l a \( l - 7 5 )) ^ ) 
(a 

* W^KT) ( -' nn < V) W 

x (in') uf(k-) . (21) 



The propagator of the off-shell r will give rise to a 
dependence of the amplitude on the off-shellness A 2 = 
El_- V l-Ml_. 

In the limit of infinite mass of the IF-boson and van- 
ishing electron mass, the unpolarizcd squared amplitude 
for the process depicted in Fig. [TO] is [24] 

I A,f|2 77 "QED "weak 

x Pvi'P~ (nJr M T - k + -k~ + k + -p + k~ -p~ + k + -p~ k~ 




The (differential) cross section for this process is pro- 
portional to the above squared amplitude, integrated over 
the threemomenta of the outgoing r-neutrino and lepton- 
antineutrino. This is a standard computation that can be 
incorporated in the computerr code if needed, yet it is 
clear that the resulting cross section will depend on the 
off-shellness A 2 and the virtual r mass-energy dispersion 
relation. In particular, at large off-shellness, the cross sec- 
tion behaves as 1/A 2 and the constant multiplying this 
parametric behavior can be fit to data. 
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Abstract. We study the multiple solutions of the truncated propagator Dyson-Schwinger equation for a 
simple fermion theory with Yukawa coupling to a scalar field. Upon increasing the coupling constant g, 
other parameters being hxed, more than one non-perturbative solution breaking chiral symmetry becomes 
possible and we find these numerically. These "recurrences" appear as a mechanism to generate different 
fermion generations as quanta of the same fundamental field in an interacting field theory, without assuming 
any composite structure. The number of recurrences or flavors is reduced to a question about the value of 
the Yukawa coupling, and has no special profound significance in the Standard Model. The resulting mass 
function can have one or more nodes and the measurement that potentially detects them can be thought 
of as a collider-based test of the virtual dispersion relation E = ^/p 2 + M(p 2 ) 2 for the charged lepton 
member of each family. This requires three independent measurements of the charged lepton's energy, 
three-momentum and off-shellness. We illustrate how this can be achieved for the (more difficult) case of 
the tau lepton. 

PACS. 12.15.Ff Quark and lepton masses and mixing - ff.30.Rd Chiral symmetries 



1 Introduction 

Why Nature has laid down exactly three fermion families 
in the same representation of the Standard Model's gauge 
groups with universal couplings to the gauge bosons re- 
mains to be explained. Many ideas have been presented 
in the literature and we quote some to exemplify the far- 
reaching implications of any experimental progress in the 
first question of what has been called the "Fermion prob- 
lem". 

"Democratic approaches" [1] are in general based on 
the idea of equal fermion-to-Higgs Yukawa couplings. One 
then needs mechanisms to generate specific lepton or quark 
mass patterns, see for example [2]. Let us also recall the 
classic work of Froggat and Nielsen [3] that introduces the 
concept of "horizontal flavor" symmetry (nowadays often 
called Interfamily symmetry) by which the three families 
are degenerate or quasi-degenerate at some high scale be- 
yond present experimental reach, due to unknown symme- 
try. At lower energies dynamical (renormalization group) 
effects amplify any small symmetry breaking term causing 
the ratios between masses that we see in current acceler- 
ators. In essence, the idea is that each of the elementary 
particle types, for example the charged leptons (e, fi, r) 
provide a representation of the posited interfamily sym- 
metry. Since the symmetry is broken at the TeV scale 
and below, revealing it requires theoretical extrapolation 
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from current experimental data. The same authors have 
followed-up with ideas based on Anti-grand unification [4] . 

An alternative line of research assumes that the known 
elementary particles are really composite objects of more 
fundamental "preons" [5]. To date no experiment has re- 
vealed composite structure beyond the usual layers of quan- 
tum pairs involved in radiative corrections around seem- 
ingly point-like sources. One can argue that a new strong 
interaction binds preons and this might conceivably be 
revealed in future experimental efforts. 

Further, dynamical breaking of chiral symmetry as a 
genuine quantum effect is well known in strong interaction 
physics [6,7] and has also been invoked in flavor physics 
as an alternative to the tree-level Higgs mechanism, for 
example by Brauner and Hosek [8]. 

However no one seems to have paid attention to the 
fact that the Dyson-Schwinger equations do provide one 
with several solutions depending on the strength of the 
coupling. In this work we call attention to this point and 
its potential relevance for the fermion problem. Each of 
the three fermions corresponds in this hypothesis to a 
quantum over a different vacuum (but in the covariant 
formalism one pursues the study of correlation functions, 
here propagators, and sidesteps the issue of the vacuum 
wave functional). Far from attempting a complete theory, 
we will show the general features of the mechanism within 
a simple model of a fermion field coupled to a scalar boson 
through a Yukawa coupling. 
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From this prospective, different flavored fermions, one 
corresponding to each family, are quanta of the same field, 
so the lagrangian can be more economically written. Each 
fermion is an elementary one-particle excitation on top 
of a vacuum that is a local minimum of the Hamiltonian. 
The excitations over the ground-state vacuum provide the 
lightest fermion family. Finally, at fixed coupling g, the 
Dyson-Schwinger equation has a finite number of solu- 
tions, so a finite number of families arise. The question on 
why three families becomes only a question of the value 
of g, the interaction coupling, and may have no special 
significance. 

The equivalent of the propagator Dyson-Schwinger equa- 
tion in a non-covariant framework is the well studied mass 
gap equation of potential models of QCD. The equation 
has been numerically solved for the harmonic oscillator 
potential [9], [10], the linear potential [11], and the lin- 
ear plus Coulomb potential [12,13]. Bicudo, Ribeiro and 
Nefediev have systematically studied the excited solutions 
of this equation [15]. Our results seem compatible with 
this prior work. 



2 Dyson-Schwinger equation and its 
numerical solution 

We examine a simple Yukawa theory for one fermion field 
coupled to a real scalar field with Lagrangian density in 
Euclidean space 

C = &(i i> - m 9 )V + g$*<j> + - m\)4> . (1) 

The Dyson-Schwinger equation for the fermion propa- 
gator in this theory is represented in figure 1 in the rain- 
bow approximation (neglecting the vertex dressing). 

To examine generic features of dynamical chiral sym- 
metry breaking we can further ignore the running of the 
wavefunction renormalization and set Z(p 2 ) — 1, a con- 
stant. This leaves one scalar equation for the fermion mass 
function that reads 

r°° M(a 2 ) 
M(p 2 ) = m*Z 2 + J o dqF( q ,p) q2 + ^ 2 , iq2) . (2) 

In the chiral limit m& — > this equation becomes homoge- 
neous and accepts a chiral-symmetry preserving solution 
M = that continuously deforms into a soft-running mass 
form when is not zero. But for strong enough kernels F 
it also admits other solutions that break chiral symmetry. 
The kernel in this simple Yukawa theory is 

F(q,p) = ^ f V^dy- I— - (3) 

4tt 3 J_j_ (q-pr+m^ 

with y the cosine of the polar angle in Euclidean four- 
dimensional space from p to q. 

To obtain excited solutions for M as function of p 2 
we apply an iterative linear method with different initial 
guesses with or without nodes for M{p 2 ). This method 



proceeds by examining linear deviations of the exact so- 
lution 

M(p 2 )=M (p 2 ) + e(p 2 ) (4) 

that yield 

poo nj / 2\ 

= m*Z 2 M (p 2 ) + ^ dqF( q ,p) q2+ ^> q2) . 

This can in turn be discretized (simply discretizing p) and 
cast as a linear system 

A l0 e = [RHS)i (6) 

and solved with standard numerical linear algebra tools. 
The kernels require only the calculation of 2-dimensional 
integrals that are also standard fare in modern computers. 

By naive dimensional analysis one sees that the inte- 
gral needs UV regularization. We accomplish this by a 
simple cut-off method. Then the parameters should be 
chosen to run with the cutoff g(A), m^(A) to ensure that 
the mass functions turn largely independent of A. In our 
results we display the dependence of the solution with g 
at fixed A. The boson mass is chosen to be 1 and sets 
the unit of the theory. 



3 Numerical results and discussion 

In this section we present the chiral-symmetry breaking 
conventional and new excited state solutions. By running 
our computer program incrementing g we meet critical g 
values go, gi, ---gn above which there are exactly n + 1 
solutions, the last with n nodes, that we find numerically. 
For example the first solution appears for a critical go — 
6.8 (of course this number runs with the cutoff). 

We show in figure 2 the three solutions obtained with 
g = 14, A = 200. 

Then in figure 3 we plot the nodeless solution at fixed 
cutoff incrementing g sequentially to show the dependence 
of the mass function on the coupling. 

Only the nodeless solution has been widely used in 
past literature for its applications in hadron physics. The 
others are sometimes rejected [16] on the basis of "wrong" 
UV asymptotic behavior, arguing that they do not match 
with conventional mass running in perturbation theory, 
and artificially imposing a boundary condition that dis- 
cards all but one solution. 

We point out here that this mismatch poses no prob- 
lem. One should interpret these solutions as the one-particle 
excitations above different extrema of the Hamiltonian or 
vacuum replicae. Therefore they cannot all be matched to 
the same perturbation theory around one particular vac- 
uum. 

If one insists on viewing them from the one vacuum 
that connects smoothly to perturbation theory, then they 
have to be written down as complicated collective states 
in terms of the one-particle solutions over this vacuum. 
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-1 -1 

• = > + ' * 1 

Fig. 1. Dyson-Schwinger equation for the fermion propagator in Yukawa theory, in the rainbow approximation. 

Comparison of the three solutions 

g=14,m=l,A=200 




50 100 150 200 

p (arbitrary units) 



Fig. 2. We plot together the three solutions to the Dyson-Schwinger equation for a fermion-scalar Yukawa coupling g — 14, 
cutoff A — 200, and the scale set by the scalar mass being fixed to = 1. Note the equation for M in the chiral limit (zero 
current fermion mass) cannot distinguish M from —M so the global sign of any one solution in this graph is irrelevant and for 
display purposes only. In particular all M(0) can be taken as positive. 

In the particular case of the Yukawa theory, the an- 
gular integral in eq. 3 can be performed analytically [14], 
but we prefer keeping a two-dimensional integral thinking 
of future work with more complicated model Lagrangians. 

In agreement with the findings of [15], we see that the 
higher excited vacua in this covariant Euclidean formu- 
lation present zeroes in the mass function. A version of 
the Sturm-Liouville theorem must be at work (this is not 
surprising since the Hamiltonian is hermitian). 

Also note that the values of M(0) generated in this 
covariant Yukawa theory show quite some hierarchy. In 
figure 2 one can see that they are well-spaced. Therefore 



this field-theory mechanism might conceivably be at work 
in the fermion flavor problem. However M(0) seems to be 
larger for the nodeless solution than for the excited ones 
(this puzzling property also appears in the equal-time ap- 
proach of [15], and the extent of its model dependence 
needs to be further investigated). To confirm this point, we 
choose a different set of parameters and implementation 
of the computer code. With — 100, a cutoff A = 10 4 
implemented as a Gaussian fall-off of the integration mea- 
sure (as opposed to terminating the grid), and g = 50, 
with a small bare fermion mass = 5 • 10~ 4 , one ob- 
tains the three solutions depicted in figure 4. This graph 
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Fig. 3. Nodeless solution for increasing values of the coupling constant (m^> = 1 and A — 200 are held fixed). Note how M(0) 
grows steadily with g (the sign of M is irrelevant). 



confirms that the nodeless solution is higher in mass, and 
that it is possible to generate sizeable mass hierarchies 
between the solutions. We use this run to also illustrate 
the dependence of Ma with the cutoff, displayed in figure 
5. As can be seen, the excited solution does not decrease 
monotonically as g — » 0, it rather seems to become a sin- 
gular solution at the critical coupling. This technicality is 
however probably irrelevant for physical applications. 

In hadron physics the strength of the interaction being 
fixed by extensive phenomenology, there would be no free- 
dom to vary g. In the physically relevant cases studied in 
[15], excited solutions exist at physical values of the cou- 
pling. In [17] Ribeiro and Nefediev have further extended 
their original results to include bound states (mesons) con- 
structed with the replicated quasiparticles. In particular 
there are Goldstone bosons over the replicated vacuum, 
that are of course not massless pions but would appear as 
excited pion states (separated from the ground state pion 
by the mass gap between the two vacua). Whether any 
conventional mesons accept a more convenient description 



as replicated mesons instead of excited mesons over the 
standard vacuum remains an open question. 

Working in Euclidean space as we do provides fastly 
convergent integrals. This technicality may be avoided by 
employing a Lehmann representation [18] that allows a 
direct solution in Minkowski space with similar results. 

4 Experimental signature 

In a free field theory, a field quantum of mass m satisfies 
the equation 

E = VP 2 + m 2 (7) 
that causes a pole in the free propagator 



In an interacting field theory this pole will be shifted 
due to renormalization from its bare to its physical po- 
sition. The dependence on the renormalization scale can 
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Multiple gap functions solving the same DSE 

(shifted for clarity) 
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Fig. 4. Three solutions to the gap equation in Yukawa theory obtained with = 100, A = 10 4 , = 5 ■ 10~ 4 , g = 50, 
confirming that the nodeless solution is higher in mass, and that a sizeable hierarchy of masses is possible. To represent them 
in a log-log plot, Mi and M2 being negative in part of their domain, the solutions have been shifted by a constant amount to 
positive values. 



be traded by a dependence on the particle Euclidean 4- 
momcntum so that 

p 2 - M(p 2 ) 2 K ' 

yielding the transcendental equation for the pole position 

E = ^p 2 + M(E 2 -p 2 ) 2 . (10) 

Another way of visualizing the function M(p 2 ) of the 
interacting theory is to trade p 2 by the space-like part 
p 2 as argument of M, generating therefore a non-trivial 
dispersion relation for virtual particles, and the zeroes of 
the mass function M appear as points where E = |p| that 
would otherwise only be reached asymptotically at large p 
(see figure 7). Tests for a non-trivial (real particle) disper- 
sion relation have been proposed and some carried out in 
the search for violations of Lorentz invariance [19]. How- 
ever these are very indirect and usually performed at low 



energy, so an accelerator-based test of the (virtual parti- 
cle) dispersion relation by separately measuring E, p and 
the off-shellness is preferable. In our case, real particles 
still provide representations of the Poincare group and 
have constant, physical mass. However virtual particles 
off their mass shell will display the running mass function 
and this can be captured by analyzing the amplitudes of 
physical processes in perturbation theory. 

The corresponding dispersion relations are plotted in 
figures 6 and 7. 

Lacking a real theory there are two possibilities we can- 
not discern at present. Either all three families of fermions 
arise as solutions spontaneously breaking chiral symme- 
try, in which case two of them should already have a zero 
in their mass function, or the first generation reflects the 
trivial solution only slightly modifying the bare mass, and 
the second and third generation solutions with dynam- 
ically broken chiral symmetry, in which case only one of 
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Dependence of M with the coupling g 
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Fig. 5. Dependence of M2 with the coupling constant at fixed m^> = 100, A = 10 4 , = 5 • 10 -4 . Upon comparing with figure 
3 we see that the qualitative behavior when the solution approaches the critical coupling and is about to cease existing, here 
g = 43, is quite different for the ground and excited solutions. 



them would present a node. In either case a zero is present 
in the mass function of one of the fcrmions. 

An interesting experiment is therefore a measurement 
of p, E and the off-shellness A 2 = p 2 — M 2 (p 2 ) of a r, 
at various momenta p. The B factories have accumulated 
large r samples [23] that are under analysis and can be 
employed for this purpose. 

The process one may investigate is depicted in Fig. 10 
in the appendix below. An electron-positron pair collides 
and annihilates into a t~ t + pair, where the r~ is off its 
mass shell. Its four-momentum can be infered from the 
center-of- momentum energy of the incoming e~e + pair 
and by reconstructing the energy and momentum of the 
on-shell r + from its decay products, as detailed shortly. Fi- 
nally, one needs the off-shellness, but this can be obtained 
from the number of counts with given E, |p|. This is ob- 
vious since the cross section for a specific process where 
the off-shell t~ decays into for example l~Div r , depends 
on the off-shellness A 2 of the r~ (see also App. A). 



Let us now detail the possible analysis: 

In the center of mass frame of an e~e + collision one 
can identify two back-to-back hadron jets, tagging the 
flavor by demanding that one of them kinematically 
reconstructs an on-shell r. This is rendered difficult 
by the undetected neutrino, that we assume to be the 
only unreconstructed track on the left side of the event. 
Complete reconstruction of the energy and momentum 
on the left side of the event (see Fig. 8) is possible with 
a vertex detector (giving the direction of motion of 
the on-shell r). The total energy is taken from collider 
calibration, and matched to the energy of the visible 
tracks. Balancing energy provides the energy of the 
missing neutrino, and automatically the momentum of 
the left side, that can be tested against the hypothesis 
of a physical r being produced. 

Next one examines the right side of the reaction, where 
the secondary and primary vertex seem to coincide (as 
appropriate for a virtual, off-shell r decaying rapidly). 
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Dispersion relation for a fermion with running mass 

Arbitrary units 



E(p)=(p +m ) 

E(p)=(p 2 +M(p) 2 ) 1/2 

E(p)=p (free massless theory) 




50 100 150 200 



P 



Fig. 6. The virtual particle dispersion relation at fixed cutoff scale for the nodeless M(p 2 ) solution. Note that the interacting 
theory falls below the usual Yukawa dispersion relation because the mass function drops with p. This is not very different from 
standard perturbative running with p except the effect is much larger in the case of dynamical chiral symmetry breaking. 



The total three-momentum is known a priori. One bal- 
ances momentum to obtain the momentum taken by 
the undetected v T , and obtains the total energy. Now, 
it is not necessarily true that E 2 — p 2 = m 2 . 
4. The number of counts with given E, p on the right 
side is normalized to a cross section in terms of the 
off-shellness <i(A 2 ) and used to obtain M from the 
value of the intermediate propagator in perturbation 
theory. 

An alternative possibility is to carry out the measure- 
ment on the expected 0(20-30) r events at the OPERA 
experiment [21] in Gran Sasso National Laboratory. Given 
the momentum of the neutrinos from the CNGS beam 
one can attempt full kinematical reconstruction of (E, p). 
However here the difficulty resides in the off-shellness since 
a small value will be forced for the identification of the r. 

The experimental test for the muon and the electron is 
even simpler, for example through Compton and reverse 
Compton scattering where the fermion is off-shell in the 



intermediate state. For quarks one needs to take into ac- 
count that they always appear in bound states, and there- 
fore the mass function is always under an integral sign and 
reconstruction from experimental data difficult. However 
a particularly simple case is 

e~e + — > 7* — > ss 

where the strangeness be completely tagged by counting 
all hyperons in both jets. 

In principle the same analysis carried out here could 
also be undertaken for Majorana particles but this is out 
of our scope in this preliminary work. 

5 Conclusions 

We have called attention to an interesting feature of held 
theory, namely the possibility of having several solutions 
of the one-particle Dyson-Schwinger equations for a broad 
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Dispersion relation for a fermion with running mass 

Arbitrary units 




P 



Fig. 7. The virtual dispersion relation at fixed renormalization scale for the solution with two nodes (we find similar result 
for one node). As the running mass actually has a zero, the dispersive curve touches the massless limit that would otherwise be 
reached only asymptotically. 




Fig. 8. In center of mass frame production where p r + = one can fully reconstruct the decay of both r leptons if only the 
neutrinos escape undetected. The analysis steps are described in the text. 
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class of theories. The excited solutions appear upon in- 
creasing the coupling constant sufficiently, breaking chiral 
symmetry. 

The spectrum of states on top of each of the vacua 
is called "replica" or "recurrence" and provides a mech- 
anism that might be interesting for the theory of flavor. 
The number of recurrences arises from no special symme- 
try, since non-linear equations can perfectly have a finite 
number of solutions (again depending on the value of the 
coupling, this can be chosen to be three). 

These solutions have been shown in the past in a Hamil- 
tonian framework [15] and we have minimally extended 
the results of these authors to show that the covariant for- 
mulation allows similar phenomena. The excited solutions 
have mass functions that present zeroes and this reflects 
in the virtual dispersion relation (feature of an interacting 
field theory that should not be confused with breaking of 
Lorentz invariance). 

Whether these solutions will be relevant to the fermion 
family problem can be tested at B factories employing 
their large r samples directly produced. 

At this point it is a meaningless exercise to attempt 
to obtain the fermion masses in this scheme. Three data 
points Mj(0) should be used to fit three parameters, the 
current fermion mass (that we have here set to zero), the 
boson mass, and g at a given cutoff, so the current pre- 
dictive power is null. However we find this field theory 
mechanism still worth attention given the ease with which 
fermion coupling universality can be incorporated by hav- 
ing three recurrences of the same field. 

Testing the off-shell dispersion relation should allow 
to immediately discard the scenario so we believe we have 
presented a valid physical hypothesis. Of course, the suc- 
cess of radiative corrections in the standard model make 
one wish for a measurement at the highest possible energy, 
maybe at a future linear collider. 

This work has been performed in the framework of the 
research projects FPA 2004-02602, 2005-02327, PR27/05- 
13955-BSCH (Spain) and is part of the Masters thesis of 
Mr. Paramo Martin presented to the faculty of U. Com- 
plutense). TVC is a postdoctoral fellow for the Fund for 
Scientific Research - Flanders and acknowledges the sup- 
port of the "Programa de Investigadores Extranjeros en la 
UCM - Grupo Santander". 



A r-pair production in e + e collisions 
A.l On-shell cross section 



We look at the 



t+t process depicted in Fig. 9, 



where all ingoing and outgoing particles are on their mass 
shell. Using the conventions of Ref. [22], the amplitude for 
the depicted process is 



-iM = u v -(p ) {-ief) v T n+ (p+) 



■Wilis 



x v 



(fe+) (iery") uf(k-). (11) 



For the unpolarized cross section, one needs the squared 
amplitude averaged over initial and summed over final 
spins 



4 s 2 



u 3r(p 



all spins 



x vl + (fc+) 7 Ve- vf (k-)u e e - (AT) 7V«'+ 4+ ( fc+ ) 

x %+iP + )lv v >+ n >- u^-{p~) , (12) 

where the spinor indices are explicitly shown. For e + e~ 
annihilation, one has Q 2 = s. One can now rearrange the 
factors in Eq. (12). Making use of the relations 



u l s (p)u e (p) = (j>+m{)^, , (13a) 
v\(p)vt>(p) = it>- , (13b) 
one can write Eq. (12) as a product of traces 
1 e 4 



£imi 2 = 

xTr[(f +m r ) 7 ,(|i+-m T )7„] . (14) 
Neglecting the electron mass, this results in 

\M\ 2 = ^ (jfe- ■p-k+- P + + k~ ■ P + k+ ■ p- 

+ m 2 T k- ■ fc+) . (15) 

This can be written in terms of Mandelstam variables 

s = (k~ - k+) 2 = (p+ - p-) 2 = Q 2 , (16a) 
t = (jfe- - p-) 2 = (p+ - k+) 2 , (16b) 
u = (jfe- - p+ f = (p- - k+) 2 , (16c) 

as 

2e 4 



{t 2 



u 2 + 4 m 2 s 



2m 4 ) , (17) 



where at high energies, the mass terms become negligible. 
In the center-of-momentum (com) frame, the amplitude is 



]T|M| 2 = e - 2 [s 2 (1 + cos 2 (0 T )) 



+ 4m 2 T s (l -cos 2 (6» T ))] , (18) 



with 6 T the corn-angle between incoming electron e~ and 
outgoing t~ . 

The differential cross section in the corn-frame is then 
given by 



da_ 

dJ7 



a' 
4? 



4m 2 



x [s(l + cos 2 (0 T )) +4m 2 sin 2 (6» T )] , (19) 



with a = e 2 /47r. The total cross section is 



47ra 2 
3s 



Am 2 



1 + 



2mi 



(20) 
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Fig. 9. On-shell r-pair production in e + e collisions. The particles' momenta and spins are indicated, as well as the conventions 
for the photon propagator and the vertices. 




A. 2 Half off-shell cross section 

In this part of the Appendix, we look at the e~e+ — > t~t + 
process, where one of the T-particles is off its mass shell 
and decays into a lighter lepton I (electron or muon). We 
discuss the case for an off-shell t~ . The case for an off- 
shell t + is clearly completely analogous. The process is 
displayed in Fig. 10. 

The off-shellness of one of the r's changes the kinemat- 
ics. The total energy in the reaction is unevenly divided: 
E T + 7^ E T - ^= VFcom/2, where W com — yfs is the total 
energy in the com-frame. The amplitude for the process 
depicted in Fig. 10 is 

- iM = u l ~ (p,) (-i ^ Y\ (1 - 7 5 )) ) 



x (-i) 



m: 



, , ■ ^ (^)H^7^(l-7 5 ) 



Hey) <+V) 1 ig » v 



{t>--M T -) ■ ' " - ' V q* 

x vf.(k+) (le^) uf(k-) . (21) 



The propagator of the off-shell r will give rise to a 
dependence of the amplitude on the off-shellness A 2 = 
E 2 -pl-M 2 . 

In the limit of infinite mass of the 1/F-boson and van- 
ishing electron mass, the unpolarized squared amplitude 
for the process depicted in Fig. 10 is [24] 



ol2 4 2 2 

El ka\2 L 77 QED "weak 

x ^Pe t -p~ (m T M T - k + -k~ + k + -p + k~-p~ + k + -p~ k~ 

A 2 i 
- — ( k+ -p Pl k~ ■ p+ + k+ ■ p+ k~ -p Pl ) . (22) 



The (differential) cross section for this process is pro- 
portional to the above squared amplitude, integrated over 
the threemomenta of the outgoing r-neutrino and lepton- 
antincutrino. This is not a simple task, yet it is clear that 
the resulting cross section will depend on the off-shellness 
A 2 and the virtual r mass-energy dispersion relation. 
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